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Abstract 

We propose a route to all-graphene integrated electronic devices by exploring the influence of 
strain on the electronic structure of graphene. We show that strain can be easily tailored to generate 
electron beam collimation, ID channels, surface states and confinement, the basic elements for all- 
graphene electronics. In addition this proposal has the advantage that patterning can be made on 
substrates rather than on the graphene sheet, thereby protecting the integrity of the latter. 

PACS numbers: 81.05.Uw,85.30.Mn,73.90.+f 



Notwithstanding its atomic thickness, graphene sheets have been shown to accommodate 
a wealth of remarkable fundamental properties, and to hold sound prospects in the context 
of a new generation of electronic devices and circuitry [1]. The exciting prospect about 
graphene is that, not only can we have extremely good conductors, but also most active 
devices made out of graphene. One of the current difficulties with respect to this lies in 
that conventional electronic operations require the ability to completely pinch-off the charge 
transport on demand. Although the electric field effect is impressive in graphene [2], the 
existence of a minimum of conductivity poses a serious obstacle towards desirable on/off 
ratios. A gapped spectrum would certainly be instrumental. The presence of a gap is 
implicitly related to the problem of electron confinement, which for Dirac fermions is not 
easily achievable by conventional means (like electrostatic potential wells) [3]. Geometrical 
confinement has been achieved in graphene ribbons and dots [4, 5], but the sensitivity 
of transport to the edge profile [6], and the inherent difficulty in the fabrication of such 
microstructures with sharply defined edges remains a problem. 

The ultimate goal would be an all-graphene circuit. This could be achieved by taking a 
graphene sheet and patterning the different devices and leads by means of appropriate cuts 
that would generate leads ribbons, dots, etc.. This papercutting electronics can have serious 
limitations with respect to reliability, scalability, and is prone to damaging and inducing 
disorder in the graphene sheet [7]. Therefore, in keeping with the paper art analogy, we 
propose an alternative origami electronics [8]. 

We show here that all the characteristics of graphene ribbons and dots (viz. geometrical 
quantization, ID channels, surface modes) might be locally obtained by patterning, not 
graphene, but the substrate on which it rests. The essential aspect of our approach is 
the generation of strain in the graphene lattice capable of changing the in-plane hopping 
amplitude in an anisotropic way. This can be achieved by means of appropriate geometrical 
patterns in an homogeneous substrate (grooves, creases, steps or wells), or by means of 
an heterogeneous substrate in which different regions interact differently with the graphene 
sheet, generating different strain profiles [Fig. 1(b)]. Another design alternative consists in 
depositing graphene onto substrates with regions that can be controlably strained on demand 
[9]. Through a combination of folding and/or clamping a graphene sheet onto such substrate 
patterns, one might generate local strain profiles suitable for the applications discussed in 
detail below, while preserving a whole graphene sheet. 
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The remainder of the paper is dedicated to show how strain only can be used as a means 
of achieving: (i) direction dependent tunnehng (ii) beam colhmation, (iii) confinement, (iv) 
the spectrum of an effective ribbon, (v) ID channels, and (vi) surface modes. 

Model - Within a tight-binding formulation of the electronic motion [10], effects of in- 
plane strain can be captured, to leading order, by considering the changes in nearest-neighbor 
hopping amplitude, t. We write t{Ri,n) = t + 6t{Ri,n), and treat the space dependent 
strain-induced modulation, 6t, as a perturbation (t ^ 3eV). It is straightforward to show 
[10] that, for smooth perturbations, the low energy Hamiltonian is 
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valid near the valleys K and K' in the Brillouin zone, with cr = [a^,ay,az), Vp = 3ta/2h 
and \1/ = ['?/';^(r), ip^{r), ip^,{r), ip^i^r)]^ is a spinor containing the electron fields in each 
sublattice and valley. Within each valley it has the form H = vf(t.{p — :;;^.4.), so that 
electron dynamics is determined by a Dirac equation in the presence of a gauge field yA.. This 
field stems from the perturbation to the homogeneous hopping amplitudes and is related to 
6t{R, n) via the complex vector potential (VP) A{r) = Ax{r) — iAy{r): 
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FIG. 1: (Color online) (a) Lattice orientation considered in the text. Thicker bonds have perturbed 
hopping, (b) Artistic depiction of a substrate (S) patterned with folds (F), trenches, dots and wells 
(A), upon which rests a graphene sheet (G). 
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The fact that ^ appears in (1) with its sign reversed for the K' valley guarantees overall 
time reversal symmetry. For definiteness, with the lattice orientation shown in Fig. 1(a) we 
perturb the vertical hopping by a constant amount over a finite region of width L. The 
perturbation and the associated A.{r) are 

5t{Ri, n) = 5t 6r,,o eiYi)e{L - Yi) (3a) 
A{r) = 5te{y)e{L-y)u,. (3b) 

The gauge field A. is oriented along u^., which coincides with the direction of translational 
invariance. Using units where v-p = h = 1, and allowing for the presence of an electrostatic 
potential V{r) in the barrier region, the wave equations for the K valley can then be cast 
as 

[-id^ -dy- AM] ^""(r) =[E- Vir)] ^^(r) (4a) 
[-zd, + dy- AM] i^^ir) =[E- V{r)] ^^(r) . (4b) 

In this formulation, the problem reduces to the study of Dirac electrons in the presence of 
VP and electrostatic barriers, and is related to corresponding studies of Dirac electrons in 
the presence of magnetic barriers [11, 12, 13, 14]. The profile (3b) has also been considered 
in Ref. 15 in modelling a suspended graphene sheet. 

Tunneling - We begin by analyzing the tunneling characteristics across the barrier-like 
perturbation of Eq. (3). Without compromising generality [16], in the remainder of the 
paper we shall be concerned with the situation 6t > and E > 0. We parametrize the 
wavefunction in the three regions of Fig. 1(a) as 



X^I _ ^ikxX+ikyy I j _j_ ^ ^ik^x-ikyy ' ^ 



e~ 



i4> 



^11 = (J^l "'■ j ^ikxx+iqyy _^ (J^f ^ \ ^ikxX-igyy 



■qjIII _ rp ^ikxX+ikyy 



where we have the kinetic momenta = Ecos(f), — 6t = Ecos(p, Qy = Esimp and the 
energy E'^ = k^ + ky = {kx — StY + g^. Substitution in the wave equation (4) leads to the 
following transmission amplitude [13, 15]: 

^ Q- ^^y I- sin (f) sin If 

cos{qyL) sin (p sin cp+i sin{qyL) (cos cos (/? — 1) 
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FIG. 2: (Color online) Transmission obtained from (5) for L = 10, 300, 5t = 0.1, E = 0.2 with 
(bottom) and without (top) gate potential V . On the right we depict the phase space for the 
cases where only the strain induced VP is present {Ax 7^ 0, F = 0), for a pure gate potential 
{Ax = 0, y / 0) and for a combination of both {Ax 7^ 0, F 7^ 0). The shaded (green) sectors 
represent the range of incident angles (j) which are not filtered by the barrier. 

This result is valid also for V 7^ with the appropriate substitution E ^ E — V inside the 
barrier. As pointed out in the case of a real magnetic field [17], conservation of kx requires 
Ecos(f) = 6t + Ecos(p, leading to strong suppression of tunneling for (p > arccos(— 1 + 5t/-E). 
When such condition is in effect, the internal angle has to be analytically continued to the 
imaginary axis — > zyj , — > iqy causing an exponential suppression of |Tp. Moreover, if 
6t/E > 2 tunneling is completely suppressed. These effects are illustrated in Fig. 2(a) where 
we plot |Tp for different values of L and V. 

Beam Collimation - Several aspects can be immediately identified from Fig. 2. Klein 
tunneling [3] is absent for a pure VP barrier, and if 7^ Klein tunneling might persist albeit 
at 7^ tt/2. As hinted above, for wide barriers tunneling is highly suppressed for certain 
ranges of that depend on 6t, E and V. This filtering effect for certain incidence angles is 
best appreciated by inspecting the phase space pictures shown in Fig. 2(b): the effect of Ax 
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is to translate the Fermi surface E"^ = k'^ + ky by St along the horizontal axis. Conservation 
of energy and momentum immediately leads to a sector of allowed incident angles, as drawn. 
Analogous reasoning applies for a pure electrostatic barrier (where the Fermi surface in the 
barrier changes size) or a combination of both. In all cases this geometrical construction 
immediately yields the transmission sector. It is clear that, whereas in a purely electrostatic 
barrier the transmission sector is symmetric with respect to normal incidence (0 = 71/2), in 
a VP barrier this sector always contains either = or tt, depending on the relative sign 
of E and 6t. Combining these two cases one can generate virtually any transmition sector 
with a single barrier, as exemplified in Fig. 2. 

This has immediate applicability in electron beam collimation and lensing, and the effect 
is easily amplifiable through a series of barriers [13]. In addition, the beam refracted by 
the barrier can approach or recede from the normal interchangeably by changing the sign 
of 6t. Alternatively, since the transmission sector depends explicitly on E, a suitable ge- 
ometrical configuration of barriers can be used to filter the energy of the incoming beam. 
In addition, the suppression of tunnelling for certain angular sectors leads to the appear- 
ance of a transport gap, as shown in Ref. 15. Therefore, even though small strain does not 
lead to a bulk spectral gap, the system exhibits an effective transport gap. Within the al- 
lowed sector, transmission through wider barriers is additionally characterized by a series of 
marked resonances where |Tp = 1, and electron flow is totally unhindered. Such behavior, 
usually associated with successive internal reflections, strongly suggests the possibility of 
confinement. 

Confinement - Dirac electrons are notoriously resilient to conventional confinement 
strategies on account of the Klein paradox [3]. The fact that our VP barriers exponen- 
tially suppress electronic transmission will be further explored to confine Dirac electrons. A 
confined state inside the barrier has the form 



From the wave equation it follows that E'^ = k"^ — h? = {k^ — 5tY + k"^, k^ = Ecosh-t}, 
K, = ii^sinh-i?, k^ — 6t = Ecosip, ky = Esimp. From the fact that E"^ > 0, this type of 
solution requires \kx — 6t\ < \E\ < \kx\. This is graphically represented by region 2 in 
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FIG. 3: (Color online) (a) Dispersion of the confined solutions, En{kx) for 5t = 0.1 and L = 200. 
(b) A density plot of the transmission t{E, k^) from Eq. (5) is shown in the upper region (region 
1), which is plotted together with the same data as in panel (a), (c) Regions 1, 2 and 3 discussed 
in the text, (d) Linearly dispersing modes for the narrow barrier discussed in the text. 

Fig. 3(c). In addition, continuity of the wavefunction requires that 

cotikyL) sin ip sinh = 1 — cos ip cosh 'd (7) 
be satisfied. When solved for Eq. (7) yields a discrete spectrum of energies for each value 
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of kx'- En{kx)- Region 2 shown in Fig. 3(c) is therefore characterized by the emergence of ID 
channels (quantum wire) confined to the barrier region, and dispersing along k^- A particular 
realization of this effect is shown in Fig. 3(a). It is clear from this figure and Eqs. (6,7) that 
such states share all the features of the ID modes typical of graphene nanoribbons. In 
particular, at the displaced Dirac point the spectrum scales as En{kx = St) ^ [n + l/2)7r/L 
and, for all purposes related to the electronic states and spectrum, in region 2 this system 
behaves as a nanoribbon. This includes the energy scales associated with the confinement 
gaps. For example, the gap is roughly Eg ^ 1/L eV.nm, so that L ^ 20nm yields a gap 
of ~ 50 meV. This is valid as long as 6t ^ Eg, which should be the case for 6t/t ~ 10%. 
Hopping variations of this magnitude can be achieved with a local strain of ~ 5% [18]. 

Surface Modes - The similarity of a VP barrier with the physics of a nanoribbon achieves 
its fullest if one notices that surface states are also possible. Just like the edge modes of 
zig-zag nanoribbons, we consider a state localized at the barrier edges. Such state should 
decay to both sides of ?/ = and y = L simultaneously. Consequently we can construct 
its wavefunction from Eq. (6) through analytical continuation of ky to the imaginary axis: 
kx — > in'- The energy will be E"^ = k^ — = {kx — St)"^ — k'^, which restricts the space of 
solutions to the areas outside both Dirac cones [Fig. 3(c)]. The solution for the wavefunction 
so constructed leads to a quantization condition analogous to (7) but, since the circular 
functions are converted to hyperbolic, will admit only one solution, valid precisely when 
< kx < St [region 3 of Fig. 3(c)]. In region 3 we then have a single state, whose dispersion 
is shown in Fig. 3(a). This mode smoothly merges with the lowest state in region 2 and, 
in the limit of L oo its energy is given hy E ^ 2{6t — kx) e~^^^^~''^\ making clear 
that this mode's energy decays exponentially from the shifted Dirac point towards kx = 0. 
Inspection of the solution in this limit reveals that the amplitude in one sublattice becomes 
much smaller than in the other, in complete analogy to the edge states in a wide nanoribbon. 
It is straightforward to demonstrate that such surface states persist in the limit L = oo where 
the problem reduces to a VP step. For a step, the surface modes have strictly zero energy, 
and occupy only one sublattice [19]. 

The results shown in Fig. 3(a) are symmetric with respect to E = [16]. This figure shows 
that the region 2U3, lying outside the continuum of the bulk of the system, supports spatially 
confined ID modes, which have the same characteristics as modes in a ribbon, including the 
presence of surface states. We can now coherently interpret what happens in region 1, 



8 



in connection with the tunnehng calculations described earlier: just as in a conventional 
Schrodinger-like barrier, the scattering states in region 1 (the bulk continuum) should "feel" 
the presence of the confined solutions inside the barrier. This manifests itself through the 
transmission resonances already shown in Fig. 2(a). In fact, when the transmission (5) 
is plotted in the {E,kx) plane, we obtain the result shown in region 1 of Fig. 3(b). The 
darker regions of this density plot correspond to the transmission resonances, which lie in 
the extrapolation of the confined modes into the bulk continuum. 

Narrow Limit - The approach used so far begs revision for the case of very narrow and/or 
large perturbations in the hopping. The clearest example is the case where L approaches the 
lattice spacing, so that hopping is perturbed only in one unit cell along Uy. Clearly, if 6t = — t 
the upper and lower half planes become disconnected and no tunneling is expected. This, 
however, is not captured by Eq. (5) which assumes overall continuity of the wavefunction 
and small perturbations. Since the narrow limit is of interest, e.g. in the case of a tight 
graphene crease, we tackle it by analyzing the localized perturbation 

StiRi, n) = 5t 6r^,oSY,,o A{r) = hvp f 6{y)u^ . (8) 

One consequence of the first order nature of the Dirac equation (4) is the discontinuity of 
the wavefunction in the presence of this potential. The presence of 6{y) in (8) imposes a 
boundary condition at the origin. In our case, there is the peculiarity that the weight of 6{y) 
cannot be even for both sublattices: for the geometry of Fig. 1 (a) the appropriate integration 
of this Dirac-delta is ^ 

r dyij^'''{x,y)6{y)=ij^^^{x,0^'-) (9) 

leading to the boundary conditions iIj^{0~) = riilj^{0~^) and ip^i^O'^) = r]ip^{0~), with 
rj = 1 + 5t/t. The wavefunction is naturally discontinuous and, moreover, for 5t = —t 
{rj = 0) vanishes in different sublattices at the upper and lower half-plane. Notice that in 
the case 6t = —t the system is comprised of two independent semi-planes with a zig-zag 
termination. Our BC reproduces the correct BC for these two zig-zag semi-planes [10, 20]. 
The transmission and reflection amplitudes as a function of (p = aj:g{kx + iki) read 

1 — 7]"^ 1 — e"^**^ 

^ ^2 _ ^-2i(t> ' = V ^2 _ g-2j0 ' (-'-^) 

with |i?p-|-|Tp = 1, and T = for rj = as expected. The interesting fact about such narrow 
barriers is that they still support surface modes. Applying the same procedure outlined above 



to the case (8) one straightforwardly obtains a state that decays exponentially to both sides 
of the axis y = 0. Its dispersion is given by = ±fs|fca;| where = 2|?7|/(l+?7^) and such 

states exist for > {k^ < 0) when rj^ > 1 (jf < 1). This is a very interesting situation: 
the perturbation of a single row of hoppings leads to the emergence of linearly-dispersing 
ID modes that live along the perturbed line. These modes detach from the continuum for 
small 5t {E < ±|/cj,.|), and reach zero energy at 6t = —t, just as expected for the two zig-zag 
semiplanes that result in that limit. This effect is illustrated in Fig. 3(d). 

In summary, our approach demonstrates that strain-induced gauge fields can be tailored 
to generate confined states, quantum wires and collimation in graphene. These results, 
together with the fact that strain has been reliably controlled in graphene [9], opens an 
exciting prospect towards all-graphene electronics. 
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